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Synopsis 


Structural vibration of a car body, caused by the engine propagates as structure borne 
sound from its source to the hood. The problem lies in controlling the sound radiation by 
controlling the vibration of those modes which couple effectively to the acoustic media. 

An experimental investigation has thus been performed to determine the feasibility of 
implementing pclyvinylidene flouride (PVDF) as error sensors on the surface of a car hood in 
a Modal filtered feedback approach to minimise the acoustic radiation . To this end, three 
PVDF sensors are positioned on the car hood such that the dominant observed responses are 
due to those modes which are more efficient acoustic radiators. Distributed PVDF actuators 
are utilised to control vibration of these modes. As a pre-test experiment, the Modal filtered 
feedback algorithm is first implemented to control the vibration of a simply supported beam 
with four actuators and sensors. The actuator and sensor placements are decided based on the 
placement and performance indices which are in-effect dependent on the system and modal 
norms, outlined by Gawronski in [4]. The modal controller is implemented on a digital 
controller board by downloading the algorithm from Matlab simulink. 

Further more, an active-passive hybrid treatment for control of vibration and sound 
radiation is investigated by applying a viscoelastic passive damping coat on to the car hood . 
This approach shows much reduction is sound radiation as the passive materials are good in 
controlling modes in the higher frequency range, while active control could be utilised to 
control the modes in the lower frequency band. 

Results from these experiments indicate that the PVDF sensors and actuators show 
much promise for controlling acoustic radiation from the hood, to a large degree overcoming 
the need for error microphones in the far field. In addition, the Modal filtered feedback 
proves to be a very productive constrained control algorithm. 
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Chapter 1 

Introduction 


The thrust of current research has been to develop intelligent structures that can adapt to 
their environment such that undesired vibration or acoustic radiation can be minimized. The 
chapter focuses on the development of this field and its applications briefly. It also gives a 
general outlay of the control of distributed continuum especially focussing on the Modal 
control. 

1.1 Motivation for Adaptive Structures 

Adaptation is typically associated with the modification, change, or adjustment that 
might occur over some period of time. In biology, adaptation is considered to be a 
characteristic of an organism that makes it better able to live in its environment. The concept 
of an adaptive structure is thus biologically inspired and is motivated by the need to modify a 
structure’s response to various stimuli (disturbances) and environmental conditions (design 
uncertainties) to meet some predetermined response characteristics. Control system design 
forms a critical path in the development of an adaptive structure and essentially classifies the 
level of adaptation possible. 

The evolution of the adaptive structuies technology, and its current knowledge base, 
was made possible by the explosive growth in the microprocessor industry over the past 
decade. Digital parameters are more readily adjusted than their analog equivalents, which 
leads to a real capability to provide adaptation in an online environment. The other critical 
element for adaptive structure design has been the development of material technologies, 
primarily electromechanical materials, also well known as smart materials. Table 1.1 gives 
the names of some smart materials and their applications. 
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No. 

Smart material. 

Application. 

Example. 

1 . 

Piezoelectric materials 

As actuators and sensors 

PZT, PVDF 

2. 

Magnetostrictive materials 

As actuator 

TERFENOL-D 

3. 

Shape Memory Alloys (SMA) 

As actuator and sensor 

NITINOL 

.4. 

Optical fibers 

As sensor 

Optical fiber 


Table 1.1 : Smart materials and examples of their applications 
Piezoelectric materials have a recoverable strain of about 0. 1 % under electric field. The 
Magnetostrictive materials have a recoverable strain of about 0.15 % under magnetic field; 
the maximum response is obtained when the material is subjected to compressive loads. 
SMAs allow one to recover up to 5 % strain from phase change induced by temperature. 

1.2 Adaptive Structure Terminology 

The various descriptive terms for adaptive structures were outlined by Wada et al. 
(1990). fig 1 . 1 illustrates various types of structural entities defined as a subset of other fields 
of research. As seen in the figure. Intelligent structure (IS) was reserved for those structures 
which incorporate actuators and sensors that are highly integrated into the structure and have 
structural functionality, as well as highly integrated control logic, signal conditioning and 
power amplification electronics. This would be the final stage of the ongoing research in this 
field. 



ADS = Adaptive structure (structures having actuators) 

SS = Sensory structure (structures having sensors) 

CS = Controlled structure (ADS with closed loop control) 

ACS = Acdve structures (Structures serving some function in 

the load bearing properties of the system) 

Fig 1.1 : Terminology of adaptive structures 
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The definition of CS is more suited for the experimental works carried out on the beam and 
hood. The active control is achieved by using piezoelectric polymers namely PVDF and PZT 
as actuators and sensors. A brief note of the piezoelectric phenomena is thus presented in 
next section. 

1 .3 Piezoelectric Phenomena 

The discovery of piezoelectric phenomena in 1880 led science and technology into a 
new dimension. The first observation of piezoelectric phenomena was done by the Curie 
brothers. Piezoelectricity is an electromechanical phenomena which couples the elastic 
(dynamic coupling) and electric (static coupling) fields. In general, a piezoelectric material 
responds to mechanical forces/pressures and generates an electric charge/voltage. This 
phenomenon is called direct piezoelectric effect. Conversely, an electric charge/field applied 
to the material induces mechanical stresses or strains, and this phenomenon is called the 
coverse piezoelectric effect. In active piezoelectric structures, the direct effect is used for 
structural measurements and converse effects for active vibration controls of the continua. 
Generic theories of structures with distributed piezoelectrics have been explored by numerous 
researchers and most of the necessary theory has been outlined by Tzou in [5]. Some of the 
elastic, piezoelectric and dielectric relations have been summed up in Appendix 1 . One draw 
back of the piezoelectric materials that is worth noting is that their performance is restricted 
by breakdown voltages, hysteresis effects and limited strain rates. 

1.4 Distributed Sensing and Vibration Control 

Structures are generally distributed in nature, i.e.., structural behaviors are functions 
of time and space, hence called distributed parameter systems. Structural identification and 
control needs a close coordination of sensors and actuators. Conventional transducers and 
sensors are generally discrete, ie... they measure spatially discrete locations of continuum. A 
severe problem can thus occur when these discrete transducers are placed at modal nodes and 
on nodal lines. The same problem arises when discrete actuators are used to control the 
continua. Eventually, these modes are neither observable nor controllable by these 
transducers and actuators. Thus, in order to observe and control the continua, distributed 
sensors and actuators are highly desirable. Thus distributed piezoelectric layers can be used 
as sensors and actuators for structural control applications. 

Modal control is an approach to control a structure by controlling its modes. Modal 
control methods have come to the limelight in the last decade primarily because of the 
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availability of numerous sensors and actuators which render the control spillover ie ..coupling 
to truncated modes and unmodelled dynamics manageable. This control is popular because it 
appeals to our simplistic and intuitive notions of single degree of freedom systems. 
Mathematically, modal control derives its power from consideration of matrix de-coupling 
transformations. The Independent modal space control (IMSC) as outlined by Meirovitch [7] 
is used in the experiments with both the beam and car hood. This method suffers from a draw 
back that its implementation by point actuators implies control of reduced number of modes 
ie.. modal truncation. Another approach to control the structure is by using direct output 
feedback control, whereby sensors are collocated with actuators. Here the difficulty is in 
deciding the suitable control gains. A good sum up of the application of active modal control 
for Smart structures has been done by Inman in [18], The ability to locate Piezoelectric 
sensors and actuators almost everywhere greatly improves a system’s controllability and 
observability properties and hence allows modal control to become an useful approach. The 
modal control can be cast either in state space form or physical space form ie.. in terms of the 
physical modes of the mechanical system. Because of the difficulties in defining such 
physical modes for mechanical system, the state space approach is used in the experiments. 
However, the modes derived from the state space analysis can be thought of as being similar 
to those of the equivalent undamped mechanical system. 
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Chapter 2 

The Euler - Bernoulli - Beam 


The Beam is probably the best structure to try implementing the modal control concepts 
as it serves as a good platform to apply the same concepts to complex structures This chapter 
intends to give a brief introduction to the theoretical modal analysis of a beam. It also 
includes the derivation of the actuator and sensor equations for piezoelectrics attached to a 
beam element. 

2.1 Euler Equation for a Beam 

Although it is well know, the derivation of the Euler-Bernoulli thin beam equation is 
presented here since it outlines the basic procedure used for the analysis of distributed 
piezoelectrics actuators and sensors attached to beams. 

p(x,t) 



Fig: 2.1 Beam section with stress resultants 
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yj 


The figure above shows an element of a thin beam excited in pure bending . For the 
following derivation it is assumed that the beam cross-section is symmetrical about the center 
-line and remains normal to the neutral surface. 

The longitudinal displacement u is given by 


u(x ) = -z 


3 w 

dx 


( 2 . 1 ) 


where w is the out of plane displacement of the beam and z is the transverse coordinate of the 


beam section. 

By Hook’s law in one dimension the stress in the beam is given by 


cr(z) = -zE 



( 2 . 2 ) 


The resulting moment in the beam can then be obtained by integrating the stress distribution 
over the beam cross section such that 


M x {x) 


f 


fj/2 

cr(z)z dz, 

hi 2 


(2.3) 


which then reduces to 

M x (x) = -EIp^- (2.4) 

ox~ 

where, /= bh 3 /l 2, is the moment of inertia of the beam where b is the width and h is the 
height of the beam, respectively. Considering the beam element with the resultant shear forces 
and moments acting as shown and neglecting rotary inertia of the element, the condition of 
moment equilibrium gives: 

^^-Ax = T f Ax, (2.5) 

3x 

where M x and Tyare the the internal moment and transverse shear forces acting upon the 
element. 

Applying Newton’s law of motion in the vertical direction gives 

-p(x) + ^L = pA~. (2.6) 

dx ox 


Combining (2.6) and (2.5) the following is obtained 


3 M 

dx 2 


pA 


' w 


3 r 


+ p(x,t) 
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and finally by using the above equation and (2.6) the Euler Bernoulli equation of beams can 
be obtained as 

„ T d 4 w d 2 w 

E, J7 +pA W = ~ p{x ’ t) (2 ' 7) 

where p(x,t) is the external load or forcing function (with units of force/unit length ) on the 
beam. 


The homogeneous form of the above equation can be written as 

El- 


9 4 w 


dx 4 


+ pA-r^-= 0 


dt 1 


A solution of the form given below can be assumed 

w(x,t) = 'P(jc) T(t ) = ¥(*) Im [e im } 

with \|/(x) as the modal form-function. 

By substituting (2.9) in (2.8) the homogeneous differential equation unfolds as 


d 4v F 


dx 4 


2T = 0, 


„ 4 P A 2 

with , /l = — O) ■ 

EI »• 

The solution of the above differential equation is given by 

T'(x) = A, coshir 4 - A 7 cos /k + A 3 sinh/ix + A 4 sin x 
Writing down equation (2.8) in the dimensionless form the following can be obtained 

3 4 w d 3"w 

El ^ +ffA TF^ 

One could express the fourth order differential equation in (2. 10) as 


( 2 . 8 ) 


(2-9) 


( 2 . 10 ) 


(2.n; 

( 2.12 


^4-k*F=0, (2.12 

d? 

where, rc 4 = — c o 2 = A 4 / 4 - 

EI yy 

The dimensionless form of the solution of the differential equation now becomes 

¥(£) = A, coshx£+ Aj cosxf + A 3 sinhK^+ A, sind; (2.1- 


The constants can be found out by applying the appropriate boundary conditions of the beam 
under consideration. 



Chapter 2 : The Euler - Bernoulli Beam 


Considering 

a simply supported beam. 


'F(O) = 0 : 

A, + A 2 = 0 

(2.15) 

SF( 1 ) = 0 : 

A cosh K + A 2 cos k + A 3 sinh k + A 4 sin K = 0 

(2.16) 

T / "( 0 ) = 0 : 

O 

II 

1 

(2.17) 

II 

O 

A, cosh k — A 2 cosath- A 3 sinh k — A 4 sin/r = Q 

(2.18) 


The natural frequencies of vibration are found by substituting the above value in the equation 
(2.13) which yields : 


d X 
d? 


= n A 7t 4 sin(n^), 


(2.19) 


As sinh K^Owe have that the constants A! = A 2 = A 3 = 0, and hence A 4 sin k = 0. As A 4 * 0, 

sin k = 0. The solution of which is k = mt where n e N + . 

Thus the eigen form-function can be written as T'n (q) = sin (mtq). 

Solution of the above expression gives, 

sm(n^)[n 4 7i 4 - X 4 l 4 ] = 0. (2.20) 


This ultimately leads to the solution 

( n7t ) 4 = (XI) 4 . 


( 2 . 21 ) 


Substituting for X and solving yields 




2 „ 2 

n 7Z 


l 2 1 


EL 


pA 


( 2 . 22 ; 


2.2 Modal Summation Procedure for a Beam 

One could write the distributed force in the right hand side of the Euler-Bernoulli equation as 
a harmonic function as written below 

p i (x,t)-p l lvn{e' 0 ‘), (2.23; 


and the displacement w{x,t) can also be written down in the form 

w(x,f) = w(x)Im{e' ft '}. (2.24 

Now writing the Euler-Bernoulli equation in the in the dimensionless form with the natural 
coordinate % , 
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d 4 w 4 3 2 w , 4 s 

EIyy d^ T+l PA ^t r = l p ( x,t ^ 


(2.25) 


Substituting (2.23) and (2,24) in (2.25) the following equation is obtained 


(2.26; 


The space function part of w(%,t) can now be written down as the sum of the product of the 
amplitude and the form-function as 




(2.27) 


where q „ is the modal amplitude of vibration. 

Substituting the above equation in (2.26) the modified differential equation is 


Aa 3 4 w a 

Z J ^(E I yy -~-l i pA£l 2 %)=l 4 p z . 


(2.28) 


For a free vibrating beam the homogeneous beam equation is given by 


3 4 vx/ 


(2.29) 


with co„ as the n ” natural frequencies of the beam. 

Subtracting (2.28) from (2.29) the following equation is obtained 


OO 

'Yj n pA{6) 1 n -Q 2 Y¥ n =p z 


(2.30) 


The normal modes \|/„(Q are also orthogonal functions satisfying the relations 


1 


0 for n^m 
const for n=m 


(2.3 1; 


Hence multiplying (2.30) with \|f„(£,) and then integrating the equation over the limits 0 and 1 
it could be obtained that 


l 1 

i, pm? -fi?) J 4 1 © 4 = J pX® 4- 


(2.32; 
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In the above equation the constant is called the norm, and is denoted by A„ 

/ 

A,=jv©<if. 

o 

The modal amplitude of vibration can be written as 

1 

A J 

q = — 2 _ 

" pA(cv n 2 -Q 2 )A n 


(2.33 


(2.34) 


The damping effects of the beam can now be included by replacing the Young's modulus of 
the beam by the complex formulation E‘+iE“, with E‘ as the storage modulus and E“ as the 
loss modulus. The loss factor r\ is defined as the coefficient E‘/E“. 

Hence the equation (2.28) can be modified as 


EI n .a + iT 7) 


3X 

a<f 


■l 4 pAQ. 2 '¥ n =l 4 p 2 


(2.35) 


Including also the effect of damping into the eigen-frequency yields 

=0) n ‘'(l + i7 1) 


and therefore (2.34) can be written as 

i 


Jp, ¥.(£)# 

0 

pA (0) n 2 (l + iJ])-Q. 2 )A n 


(2.36) 


(2.37) 


Multiplying the numerator and denominator of the above equation by its complex conjugate 
one could obtain 




(6j n 2 -a 2 )-i(^ n X 

pAA[{co n 1 -Q. 2 ) 2 +co n A r ] 2 } 




0 


(2.38) 


Taking the modulus of the complex modal amplitude yields 

^(ty n 2 -ft 2 ) 2 +<y„ 4 77 




pAA(a ) n 2 -Q. 2 ) 2 +o). 


4 2 1 

njf_ f 

ajrf] J 

n / J o 




(2.39) 


2.3 Modelling of Beams with Piezoelectric Actuators 

The common form of arrangement of the Piezoelectric actuator is the anti-symmetric 
wafer configuration. In this arrangement two identical piezoelectric elements are located 
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voltage is applied across the electrodes (in the direction of polarization) the actuator induces 
surface strains in the beam through the d$ j and dyj mode of the piezoelectric material 
response. This configuration is known as the wafer arrangement since the piezoelectric 
element is very long (in the x and y directions) compared to thickness (in the z direction) 
through which it is polarised. 

One dimensional antisymmetric wafer actuators 

Important assumptions for this configuration are : 

♦> the beam is covered by a layer of thin piezoelectric material of thickness, h a which is 
perfectly bonded to the beam and strains only in the x direction. 

♦> inertial effects of the piezoelectric element are ignored which is valid if the element is thin 
and lightweight compared to the beam system. 


When a voltage is applied across the unconstrained layer (i.e.. not attached) piezoelectric 
element the actuator will strain by an amount £ pe in the direction 1 which is parallel with the 
axis as dictated by 


£ 


P e 


h.. 


(2.40) 


where U a is the applied voltage in the direction of polarization , h a is the actuator thickness 
and d 3 i is the piezoelectric material strain constant. 

By the assuming that the strain distribution is linear as a result of Kirchoff s hypothesis of 
laminate plate theory it can be written that 


£ A (z) = Cz +£ 0 , (2.41) 

where C is the slope and e 0 is the z intercept. The above equation can be decomposed into 

the sum of a antisymmetric distribution Cz - due to the flexural component - about the center 
of the beam and a uniform strain distribution e 0 - due to the longitudinal component . It 

must be stated here that the longitudinal component is present only in the case of asymmetric 
wafer arrangement with only one actuator attached to one side of the beam. In case of 
antisymmetric wafer arrangement, which was investigated in this research, pure bending 
occurs and the longitudinal component vanishes, i.e. e 0 =0. The figure below shows the 
above stated decomposition for a asymmetric wafer configuration. 
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Flexural Longitudinal 


Fig : 2.2 Asymmetric wafer configuration 

With the linear strain distribution and Hooke’s law , the stress distribution within the beam is 
given by 

a h (z) = E h (Cz + e 0 ) (2.42) 

where E b is the Young’s modulus of the beam material. The stress distribution within the 
piezoelectric actuator <J pe (z) is a function of the unconstrained piezoelectric actuator strain, 
the Young’s modulus of the actuator material E pe . 

Thus the stress distribution can be written as 

<T P '\z)=E r ,(Cz+e „ + £,,) ■ (2.43) 

= + «„-«„) • C 2 - 44 ) 

Now applying moment equilibrium about the center of the beam, 

— h hj-, hfo +h tt 

j<j pr 2 (z)zdz + | a h (z)zdz+ zdz = 0 (2.45) 

~ h i, ~ h „ ~ h i, h i> 

where the superscripts 1 and 2 denote the piezoelectric elements as shown in the fig. 2.3. 
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z 



Fig : 2.3 Antisymmetric wafer configuration 


After substitution of (2.42), (2.43) and (2.44) in (2.45) the strain slope C is obtained as 


C = K f s 


pe 


(2.46) 


where the material-geometric constant is specified by, 


K f 


3 E pe [Qh,+K? -K] 


nEM+hj -V]+w> 


(2.47; 


The line moment m x induced in the beam by the actuator is again given by 

m x (x) = E h IK f £ pe . (2.48) 

Now writing the Euler-Bernoulli equation of the beam in the moment form including the 

actuator induced moments, the beam - actuator equation of motion is given by. 


d 2 [M X (x)-m x (x)] 


dx 2 


CD 2 pAw = 0 


(2.49) 


where M x is the internal beam bending moment as shown in fig 2.1 and m x (x) is the actuator- 


induced line moment. 


For a finite element length L a of actuator the actuator-moment equation can be modified as 

m x (x) = C 0 e P '[H(x)-H(x-L a )} (2.50] 
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where H(. ) stands for the unit Heaviside step function defined as 


H (x ) 


I , x > 0 

0 , x < 0 


and 


C 0 =£*/# 


The above equation implies that the induced moment only exists at every point under the 
location of the finite actuator. Substituting the moment equation into (2.49), taking the 
second partial derivative with respect to x and moving the actuator terms to the right hand side 
of the equation one find 

- ^ 2 (A) -co 2 pAw= C 0 £ pe [S\x) -S'ix-LJ] (2.51) 

where S‘(.) represents the derivative of the Dirac delta function with respect to its argument. 

Hence by substituting the force per unit length p z (x) in the Euler-Bemoulli equation 
corresponding to the actuator force as 

P ; W = C 0 £ pe [S' (x) -S'(x- L u )] (2.52) 


2.4 Modelling of Beams with Piezoelectric Sensors 

Point sensors are usuaily employed in arrays of transducers whose outputs are processed 
to obtain some estimate of a required variable or state of the system to be controlled. The 
basis of this approach is that in order for the control to be effective with out observation 
spillover, the controller has to be designed to observe only those modes which are to be 
controlled. The main disadvantage of this approach is due to the signal processing 
requirements necessary to process the transducer outputs and thus obtain estimates of the 
required variables. The output of the sensor is sensitive to the accuracy of the shape of the 
sensor as well as its positioning on the structure. 

Considering a one -dimensional thin beam covered with a thin layer of piezoelectric material . 

❖ Its assumed that the piezoelectric material is mechanically isotropic and that the sensor 
has constant properties along its length and is also thin compared to its length. 

❖ It is also assumed that the sensor has no effect on the motion of the beam , i.e. the stiffness 
and inertial force of the sensor are very small compared to the stiffness of the beam. 
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For a piezoelectric material the electrical and mechanical constitutive equations are coupled 
and in one-dimension they can be written as 


— ?;,+«, ( 2.53 

D,=d M T n +ejE 3 . ( 2.54 


where, 

D 3 is the charge per unit area expressed in Cm' 2 , 

£3 is the electric field applied in the poling direction expressed in Vm' 1 

e 3 3 is the permitivity under constant stress. 

^31 is the piezoelectric constant and relates the strain to the electric field E in the absence 
of mechanical stress. It is expressed in mV' 1 
Tu is the stress in the material 

s E is the compliance of the piezoelectric material when the electric field is constant 
Sii the strain in the material 

( 2.55 
( 2.56 


The above equations can be transformed into 

■Tl =< 1I 

D 1 =e iI 5 ’ ]I +£ ! /(l-/c 2 )£,. 


where 

F F 2 

c = l/s\ 1 J is the Young’s modulus under constant electric field in Nm' 

p 

£31=6(3 i/s 11 " is the constant relating the charge per unit area to the strain for short circuited 

-2 

electrodes in Cm 

k 3I 2 = d :n 2 /(s ll E £„ 1 ) is called the coupling coefficient of the piezoelectric material, a high 

value of which is desirable for efficient transduction. 

£3/ = e^ r (1 - k, 2 ) is the permitivity under constant strain. 

According to the Euler-Bernoulli assumption, the axial strain and the curvature are related by 

S„=(A„+A,>0 (2-57) 

where h b is half the thickness of the beam and h s is the thickness of the piezoelectric sensor. 

If U s is the voltage produced in the poling direction of the piezoelectric material by the 
deformation, then the electric field £3 can be written as 
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Substituting the above equation in (2.56) yields 

U - 6 31^11 ft .v 


e-,,S,,h. 


'33 


'33 


Defining the quantities: 


(2.59) 


— = £>-, = freechargeperareaon thepiezoelectic element. 

A 

— = eS =Piezoelecric polarizatbn chargeper area due to strain, whereA is theareaof electrode 
A 


In (2.59) the first term is Zero as there is no voltage applied to the material. Now after 
substituting the equation (2.58) in (2.59), and also including the aforesaid argument the 
following expression is obtained 


rr K / 1 , W 

- — 7 (K + K )( e 3i tt) • 

£33 dx~ 


(2.60) 


It is desired to develop a modal sensor , i.e. a distributed piezoelectric sensor whose output is 
only related to the motion of one particular mode of the beam. The motion of the simply 
supported beam can be written as 

w(x,t) = ?.(0 ¥.(*). (2.61) 

n = 1 


Substituting (2.61) in (2.60) its found that. 





n=l 




h 

^33 


3 2x W 

dx 2 


(2.62) 


Though the concept of point sensors has been used, in reality, the sensors have some definite 
size. Hence an averaging over the length of the sensor has to be carried out. Thus the 
modified expression of U s yields 


OO 

°.=I. 




K (h h +K) cd 2 %(x) 


n=l 


L 




dx 1 


■dx 


where l s is the length of the piezo sensor. 


(2.63) 


Now, this sensor voltage can be measured in two different ways either by using a resistor or 
by using a capacitor in parallel as shown in the figure below. 
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(a) (b) 

Fig 2.4 : Measurement of sensor voltages using (a) Resistance and (b) Capacitance 


For a capacitor with capacitance C, the measured sensor voltage is given by, 

V't 1 Vd 2v F 

U M =Cq p =C A e,, (\ +A,) Jjkj 


0 0 


For a resistor of resistance R, the measured sensor voltage is given by, 

A 


da v— > 1 

V u =R^RAe»(\ + K) Jj, T 


1 

0 l ' I 

When n sensors observes n modes the vector of sensor voltage yields 

U,=' Vq. 

Hence 


~U tl ~ 


U ,2 



_u _ 



'F u 
'F 2. 


*F 12 
V F 22 


'F 'F „ 


i ^ 

where for the i sensor and n mode. 


'F i„ 

'F 2 n 


*F 


?i 

<?2 




f «2 - l2 


'F 


-VA A) 


"33 


«> 




e " 1 3x 2 * 

V Xl > 


(2.64 


(2.65 


( 2 . 66 ) 


(2.67) 


( 2 . 68 ) 


here xij and xij represent the coordinates of the ends of sensor i. 
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2.5 Modal Amplitudes for Excitation by a Point Sinusoidal Excitation - Shaker 
and Sinusoidal Excitation by the Actuator 

From (2.39) it is clear that the only difference in computing the modal amplitudes for a shaker 
excitation and an actuator excitation lies in the integral term containing the force p z Here an 
excitation applied on a simply supported beam is considered.. 

Distinguishing the integral for both cases, this could be written down that 
1 

/, = for theshaker. (2.69) 

0 

1 

I act - J*P- HCf X F„ (£M|, for thepiezoelectic actuator. (2.70) 

o 

In the case of a shaker, the property of the Dirac- delta function can be utilised as the 
represent the point load excitation, which leads to 

i 

7 v = J<5(£-£ 0 )Fsin(n;r£)<i£ (2.71) 

0 

= Fsin(n^ 0 ). 

In the case of the piezoelectric actuator, the force vector could be written using (2.52) in 
dimensionless form as 

P m = C„ep8{l;-tt-S . (2.72) 

Where and £ 2 represent the two edges of the actuator in the natural coordinate system. 

Hence one can arrive at an expression for the integral I acl as 

1 

I, a = { -£,)-*( # -&>]sin Qaftdt (2.73) 

0 

= C 0 e pc n7C [cos(n7r£ 2 ) -cos(n^,)]. 

This expression helps to calculate the frequency response function of the actuator excitation. 
Hence it becomes a very important tool for designing the actuator locations and maximum 
control spillover reduction. The input control function p(x) could thus be tailored to achieve 
the required modal excitation in the control field. 
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Chapter 3 

Modal State Sp^ce Formulation and Modal Control 


The design methods based on the state space approach, often called Modern methods , 
start from a system description using the first order differential equations governing the state 
variables. The formalism applies equally to single input single output (SISO) and multi input 
multi output (MIMO) systems. 


In the state variable form, a linear system is described by a set of first order linear differential 
equations, 


x = Ax+Bu+Em{ ( 3 . 1 . 

y = Cx+Du+\v 2 ( 3 .2) 


with the following notations: 
x state vector, 
u input vector, 
y output vector, 
vtq system noise, 
vv 2 measurement noise, 

A system matrix, 

B input matrix, 

C output matrix, 

D feed through matrix, 

E system noise input matrix. 



Fig 3.1 : State space form of a linear system 


The feed through component D only occurs in special cases where the input quantities, e.g. 
actuator forces and moments, have a direct influence on the sensor measurement. It could be 
a result of sensor type - e.g. accelerometer or location - when collocated with the actuator, or 
as a result of the modal truncation, the so called residual mode. The system noise w { may 
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iv 


include environmental loads, modeling errors, non-linearity’s and noise in the input vector. 

The measurement noise wj includes the sensor noise and the modeling errors. 

3.1 Modal State Space Description for a Beam with External Excitation 

By substituting for the displacement as w(x,t) = Z qJt)W(x), the Euler-Bemoulli 
equation of a thin beam can be written down in the modal coordinate q n (t) as 

4„+2 C„G) n q n +CO*q n = p nz (3.3) 


where p Ilz is the modal load and is given by, 


J p z (x,t)^ n (x)dx 

2 

pA^ 2 {x)dx 


(3.4) 


/ 

The damping was assumed to be a hystherisis one, by replacing the young’s modulus £ by a 
complex formulation E‘+iE“, with E‘ as the storage modulus and E“ as the loss modulus. 
The loss factor rj is defined as 


V=K = 


ET_ 
E ' 


(3.5 


The above second order equation implies that there are two state variables. They can be 
selected as 

* i = 

x 2 =q n (3.6) 


With this choice the equation can be written down as a set of two first order equations, 



* 2 = -2^m,x 2 -co n 2 x l + Pnz , or 


■ 

X \ 

\ x iJ 


(■ ^ 


0 1 
2 


^0 f 0 ^ 


-2 Coy n ){q n 


\Kj 


pM) 


(3.7) 


where the input factor h n can be written down as 

| p z (x)%(x)dx 

h,=^— c • (3-8) 

oA\ y ¥. 2 (x)dx 
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This can be diagrammatically represented as shown in Fig 3.1. 



Fig : 3.2 System with external excitation h 
where, the vector of input h for mode n can be written down as 



In the shortened form (3.7) could be written down as 

A 

x = Ax + h p z 


2. i 


(3.9) 


(3.10) 


3.2 State Space Form with Included Modal Actuator Input 


In the same lines as the modal input excitation vector h was described, the actuator 
input matrix B for a single mode can be written down as 



where the actuator input factor b n can be expressed as 

£ n= J . (3.12). 

pA^ n 2 (x)dx 

l 

Now including the actuator input into the modal state space form without h, the modified 
equation is 


f 


u, 


\ 


( ■ "l 
X, 




f 0 

2 



n 


i -0) 




■ 2 Co). ) 


+ 


.... 0 ' 
,b in .... b mnj 


(3.13) 
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where m is the number of actuators and u is the actuator input voltage applied separately to 
each actuator. 

This can be diagrammatically represented as shown below. 



Fig 3.3 : System with actuator input B 
3.3 State Space Form with included Sensor Output 


The output vector y, which is the signal generated by the sensors due to the deformation 
experienced during bending of the beam could now be explained more in detail. 

Using (2.65) and replacing the voltage U s by the state variable y, the following could be 
obtained 


y. = cJ2j9.0to.<t). 

n=l 

r.(x)?,«). (3.14) 


If a velocity sensor is used, then the first expression of the above equation can be got by a 
time integration of the second expression. The constant C amp is the amplification factor of 
the signal conditioner. The state variable y could thus be represented as shown below 


(y') 

=c 

'Pi- 

m 0... 

**•* X n m 

... 0^ 

Uv> 

^ amp 

U- 


u/ 

- ^nj 


< 7 „ 


(3.15; 


UJ 
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whieh can be written down in short form as , y = Cx where C is the sensor output matrix. 

The above discussion could be summed up with a block diagram which shows the system 
along with the sensor output matrix . 



Fig 3.4 : System with the sensor output matrix C 


3.4 Active Beam with Sensors and Actuators 

Now combining all block diagrams of the excitation, actuator input and the sensor 
output in to one single equation, the actual state space model of the beam could be 
formulated. The block diagram in fig 3.4 shows this. 



Fig 3.5 : Block diagram of the active beam with the sensors and actuators 

Formally the beam equation could thus be written down as 

x = Ax+ Bu + hp z 
and y = Cx. 


with the matrices A, B, h and C having the values discussed in the earlier sections. 


(3.16) 

(3.17) 











Chapter 3 : Modal State Space Formulation and Modal Control 
3.5 Modal Control 

The signals from the output y(t) weighted by the elements of the gain matrix G could 
be fed back as the input for the system, where the gain matrix is connected parallel to the 
system matrix A. 

Modal control is a term used to describe a wide variety of control techniques. In 
general, modal control is an approach to control system design in which the state variable are 
chosen such that the unforced behavior of each state variable is relatively independent of the 
behaviour of the other state variables. Here a special modal control technique is used, called 
the Independent modal space control ( IMSC ) , first coined by Meirovitch [7], in which case a 
distributed feedback control is possible. This is because the gain factor for the displacement 
(Gj ) and the gain factor for the velocity (Gf) satisfy the eigenvalue problem of the system. It 

is to be stated here that this control is only possible when the number of actuators (m) is equal 
to the number of control modes (n c ) or in other words the number of actuators is equal to the 
number of sensors. In the above stated cases the closed-loop modal equations are independent 
of each other as G[ and Gj are diagonal. If n c >m then IMSC is not possible as G\ and Gi are 
no more diagonal and hence coupled control has to be resorted to. [Meirovitch] 

3.5.1 Independent Modal Space Control (IMSC) 

The block diagram of the system and the controller is shown in Fig 3.5. The system is 
represented by the system matrix A, the actuator input matrix B, and the sensor out put matrix 
C v . The index v indicates that the sensor signals are proportional to the velocity of the 

structure in motion. The matrix M is the filter matrix which transforms the desired 
displacement w into appropriate actuator voltages {u s j. Since the object of vibration control is 
to reduce the vibration amplitudes, the values of the displacement w and the actuator voltages 
{«y} are set to zero. The output vector {y*} contains the sensor signals of n PVDF point 

sensors which are located on the surface of the structure. Since the positions, the geometry 
and the material constants of the sensors are known, the output matrix can be calculated by 

using (2.64). Assuming that the output matrix C v is non-singular, the inverse matrix C v _1 

exists. The inverse matrix C v _1 becomes the so-called mode filter/ analyser as by multiplying 
it bv the outnut vector { y * } , this yields the vector of modal velocities {q v }. Therefore a 
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system of linear equations is created as shown in (3.18 ). This type of filtering state space 
quantities by using the output vector is called output feedback. 



Fig 3.5 : Modal filtered feedback (IMSC) block diagram 
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The modal state space vector (xj contains modal displacements and modal velocities. 
The modal velocities are evaluated by carrying out a time integration of the filtered modal 
velocities. This is show in the following equation: 



0 


J 4\ dt 


f qS 

x} = 

0 

+ 

1 q-,dt 


q. 


4 1 

* 

' 

J 

0 


?i 


^4n ) 


A 


<4n, 


K U J 


( 3 . 19 ) 


The controller is defined by the proportional feedback of the modal displacements and 
modal velocities. The gain matrix G is composed of two diagonal matrices Gj and Gi- Gj 
contains the gain factors to feedback the modal displacements and the second diagonal matrix 
Gj is defined by the gain factors for the feedback of the modal velocities. As a result of their 
diagonal nature the closed-loop modal equations are reduced to independent equations and 
hence the control is called IMSC or natural control. This control has very important 
implications, as the open-loop eigen functions are closed- loop eigen functions as well. 
Hence, in IMSC the control effort is directed entirely to altering the eigen values, leaving the 
eigenfunctions unaltered. As the stability of a linear system is determined by the system 
eigenvalues, with the eigenfunctions playing no role, in IMSC , no control effort is used 
unnecessarily. 


The multiplication of the gain matrix G and the state space vector x results in the control 
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Formally (3.20) can be written down in short form as 

{ r } = [ G, G 2 ] [ x ]. (3.21) 

The elements of the control vector (rj are related to the filtered modes. Therefore, the 
controller works within the modal state space and a retransformation of the control vector {rj 
into the physical space is required. The retransformation is carried out by the so-called mode 

synthesiser which is defined by the inverse modal actuator matrix B" 1 . The multiplication of 

the inverse matrix B' 1 and the control vector {r} leads to the input vector {uj. The input 
vector {u} has as many elements as actuators are available. 



(3.22) 


The whole of the above discussion can thus be summed up in state-space form as 


{«} = B-' G cr y 

(3.23) 

and hence 


x = (A+B[B' ] G C; l C]) x 

(3.24) 

y = C x. 

(3.25) 


3.6 Determination of Modal Gains- Pole Allocation 

One fact that has to be laid down is how to determine the modal gains Gj and G 2 . 
Two of the most widely used methods, both of which have been used this study are pole 
allocation and optimal control (Fuller [1], Meriovitch [7]). In this section a brief description 
of the pole allocation method is outlined. 

3.6.1 Pole allocation 

In pole allocation, the closed-loop poles are selected in advance and the gains are 
determined so as to produce these poles. In other words, this technique is based on equating 
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the characteristic equation of the system matrix for the closed- loop system to some desired 
polynomial of same order. 

The equation for the full state feedback can be written down 

x = (A-BG T )x + Hp z . (3.26) 

The characteristic equation for the closed-loop system can be expressed as 

\sI-A + BG T \ = 0. (3.27) 

When the eigenvalues (s+X\)( s+h-i) (s+ /L„J of the desired system can be specified, it is 

possible to write a corresponding desired characteristic polynomial of the same order as 
equation (3.27): 

A^=(s+Xi)( S+X 2 ) (s+ A. lt ) =0. (3.28) 

Equating (3.27) and (3.28) we can solve for the gain G so that the equality is satisfied. 

In the IMSC, the procedure is simple. Denoting the closed-loop eigenvalues associated 

with the n !h mode by -a,, + ifi n , the solution of (3.26) can be written as 

x n (t) = c n e i ~ a " + ' /} " )r . (3.29) 

Inserting (3.29) into (3.26) and separating the real and imaginary parts yields 

and 82 n=^n n =l,2,3,...n c . (3.30; 

To guarantee asymptotic stability, however , it is only necessary to impart the open-loop eigen 
values some negative real part and it is not necessary to alter the frequencies. This can be 

achieved by letting p„ = (A.) 1/2 = co„ , where G)„ is the n th eigen frequency of the open-loop 
system. 

Hence, the frequency-preserving control gains are 

Su= a n and £ 2 „ =2 ^ w ith n = l,2,3,...n c . (3.31) 

The fig (3.7) below sketches the Pole Placement method for a single degree of freedom 
oscillator. Here it is clearly shown how the real part is shifted more to the left of the 
complex plane. 
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Fig. 3.7 Single degree of freedom oscillator: (a) Open-loop and (b) Closed-loop poles 


Though the fact that poles of the closed-loop system can be located arbitrarily in the 
complex plane is remarkable, this can always be done provided the system is controllable. 
Although, it is important to state that even if the above stated condition is satisfied, pole 
placement may not be feasible because the control effort involved would be too large, or 
because the large values of the feedback gains would render the system oversensitive to noise 
or modeling errors. 

3.7 Stability 

As stated in the earlier section, the stability of a linear feedback system can be most 
readily be determined by inspection of the position of its closed loop poles. The equilibrium 
ci iJtA v — n ic in the. sense of Liapunov if, for every e > 0 there is some 5 < such that, if 
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II xq II < $, then II x II < £ for all t < to- in this statement II . II stands for a norm, measuring the 
distance to the equilibrium. States which are not stable in the sense of Liapunov are unstable. 
A linear system is stable if the real parts of all its poles (eigenvalues of [A]) are negative i.e. 
Re(Xj) < 0, or in other words, if all poles of A are to the left of the imaginary axis when 

plotted in the complex s plane. It is worth noting that any system with a pole lying exactly 
on the imaginary s axis, i.e. having zero real part, is not stable as according to the definition 
stated in the above sentences. Meirovitch [7] has described such systems as also to be stable, 
and describing the systems satisfying the definition of stability as stated above to be 
asymptotically stable. 

It is often the case in practice, however, that the closed polynomial forms of the 
transfer functions are generally not directly available, and stability must be assessed from 
input-output measurements, made on the systems before control. The Nyquist stability 
criterion provides such a method, using the open loop frequency response. The Nyquist 
criterion states, that a closed loop system is stable only if the polar plot of the open loop 
frequency response does not enclose the point (- 1,0) in the s plane. 
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' Chapter 4 

Actuator and Sensor Placements 


An important question to ask in designing a closed loop system is where to place the 
sensors and actuators on the structure. It is also important to decide the shape of them. 
Control theory provides an answer to these questions by introducing the concepts of 
controllability and observability. A structure is controllable if the installed actuators excite all 
its structural modes. It is observable if the installed sensors detect the motions of all the 
modes. Classically, according to control theory the system is if the controllability matrix - in 
this case input matrix B - is of full rank and it is observable if the observability matrix - in this 
case output matrix C - is of full rank This approach is often not enough for practical 
engineering problems where a more quantitative information is needed. In this study the 
quantitative approach is used in terms of the controllability and observability grammarians - 
in terms of H 2 and EL norms - which represent a degree of controllability and observability of 
each mode. Thus the approach used consists of the determination of the norm of each sensor / 
actuator for selected modes which is the normalized input B and output C matrix, and then 
grade them according to their participation in the system norm. This is a computational 
effective approach with clear physical interpretation. 

The problem can thus be summarized in the following way: Given a large set of sensors and 
actuators, the placement problem consists of determining the location of a smaller subset of 
sensors or actuators such that the H 2 or EL norms of the subset is as close as possible to the 

norm of the original* set. 
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4.1 Norms 

System norms serve as a measure of the system size and in this capacity they are used 
in the actuator / sensor placement procedures. For flexible structures the H 2 norm has an 
additive property. It is a root-mean-square (rms.) sum of the norms of the individual modes. 
The HU norm can also be determined from the corresponding modal norms, by choosing the 
largest one. In addition both of these norms of a structure with multiple inputs (or outputs) 
can be decomposed into the rms. sum of norms of a structure with a single input ( or output). 
These properties allow for the development of unique and efficient actuator / sensor 
placements. 

4.1.1 The H 2 Norm 

If A , B and C be the system state-space representation of a linear system, and if 
G(w)=[C](jcofl]-[A]) ~ l [B] be its transfer function and I is the unity matrix, the H 2 norm of 
the system is defined as (Gawronski [4]) 

ll[G(ty)]ll 1 2 =— f tri[G*(Q)))[G(Q))])dC0, (4.1) 

2jc J 

where * stands for the complex-conjugate-transpose of a matrix and tr stands for the trace of 
the matrix. 


Mode Norm 


If G„(w)=[ C„](j ' l [B n ] is the transfer function of the n th mode, the H 2 norm of the 


n th mode can be derived as (Gawronski [4]) 


H[GJM 


nr/*jM [cjiia 


(4.2) 


2^0^ 

where B n and C„ arc the input and output matrices of the n th mode for S inputs (actuators) and 

R outputs (sensors). 
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System Norm 

If G(co)=[C](jco[I]-[A ]) ~ [ [B] is the transfer function of a structure then the system H 2 norm 
is approximately the rms sum of the modal norms i.e. 

(4.3) 

where n c is the number of control modes. 

4.1.2 BL Norm 

The BL norm of a stable system is defined as (Gawronski [4]) 

H[G(^]IL=To«([G(aJ]), (4.4; 

where o max ([G( co)]) is the largest singular value of G(©). The BL norm of a SISO system is 
the peak of the transfer function magnitude, see also fig. 4.1. 

Mode Norm 

th 

The BL Norm of the n mode is estimated as 

l|C ,] l5 Ma, (4 .5) 

^Cn^n 

where ^„and co„ stands for the modal damping and the eigenfrequency of the n th mode, 
respectively. 

System Norm. 

Due to the almost independence of the modes, the system BL norm is the largest of the mode 
norms i.e. 

NmsTiitGjiL. (4.6; 

This property says that the largest modal peak response of a lightly damped structure 
determines the worst-case response. 

From the above definitions, ll[G]ll 2 can be interpreted as the area enclosed by the 
transfer function and ll[G]IL is the value of the largest peak of the transfer function. This is 
shown in the fig. 4. 1 for a system with three controlled modes, e.g. n c = 3. 
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Magnitude, 
mode 2 



(a) Frequency, rad/sec 


Magnitude, 

system 



(I)) Frequency, rad/sec 


Fig. 4.1 : Hz and H{inf) norms (a) of the second mode, and (b) of the system. 
There is one more important property of these norms which make them so vital for the use in 

placing actuators and sensors. This is the additive property of a set of actuators/sensors for 

each mode. 


4.1.3 Additive Property for a Set of Actuators / Sensors for each Mode 

The two norms of the n h mode of a structure with a set of S actuators is the rms sum 
of the norms of the mode with each single actuator from this set 


II G, 


n " 2 ~ 


Z'^ 


n = l,2,...ii , 




where II G nj ll 2 = 




(4.7) 
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Similarly the two norm of the mode of a structure with a set of R sensors is the rms sum 
of norms of the mode with each single actuator from this set 


R 

"GJts 2 I[G - ]II > 2 ' n=l,2,...<> , 

I JW 

where II G nk II 

2 ^ 


(4.8) 


4.2 Actuator Placement 


The actuators in our case are distributed piezoelectric sheets which are bonded to the 
surface of the structure. The efficiency of the actuators can be expressed in scalar 
performance indices. 


4.2.1 Actuator Placement Indices Matrix 

The placement index indicates the controllability of the jth actuator with respect to the 

nth mode and is defined as : 


Jits* 3 »2 

* II [G] ll 2 
d 

H[GjJ II 
ll[G]IL 


(wherej = 1,...,S) by using theH 2 norm 


by using the EL norm 


(4.9) 


(4.10) 


The placement indices are then arranged into the placement indices matrices [E ] 2 and [E] „ 



"cr 2 .. 

.. .. <7 

l \s 
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*■ “ 



- 
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•* 
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<t 2 .. .. 

L *i 

.. <J 2 ! 

l nS _ 

i 

^oo *• •* 



Each row of these matrices is related to a mode and each column corresponds to an 
actuator. Since the output matrix is included in the calculation of the modal norms the 
actuator placement indices are dependent on the sensors. Therefore the actuator performance 
is estimated with respect to the sensor performance. If the performance of the actuator is 
required without including the sensors, then the input matrix has to he replaced by the unity 
matrix. The unity matrix does not influence the actuator norm calculated by equation (4.7) 
and emulates the ideal sensor array which enables the observation of all modes. 
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4.2.2 Control Performance Indices 

The actuator placement indices are evaluated for each actuator and each mode. In 
order to determine the performance of all actuators together a control performance index is 
introduced. The calculation of the performance index starts with the summation of all S 

actuator placement indices of the n th mode in the way of root-mean square values for o a2n or 
by choosing the largest actuator placement index related to the n th mode for cr^. 


< 7 . 



(4.12) 


= 7 • (4-13) 

These values indicate the control performance of a selected number of actuators with respect 
to a selected mode and they are collected into the modal control performance vector as 


iY/~, <T.J and (4.14) 

<7,./ (4.15) 


To obtain the numerical expression for the control performance for all modes the values <J ain 
or G a '„ n are summarized by 


oc.. 


1 ”<■ 


<J„ 


(4.16) 
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(4.17) 
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(4.18) 


B = "rj 
r a 2 nl 


c n ~\ 


n 




a 2 n ’ 


(4.19) 


n=l 


a. and a are the so - called arithmetic mean (am) control performance indices of an actuator array anc 


4 2 


and J3 a are the geometric mean (gm) control performance indices. 


4.2.3 Design procedure 

The design procedure its done by the following steps: 

1 . The modes to be controlled are selected. 

2. The actuator array is created according to the shape of the lowest mode. 

3. The control performance indices p a2 and (3 a «, are calculated. 

4. If p and p a „are equal to zero, then an actuator array is created according to the shape of 
the next higher mode. 
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5. The control performance indices cc Q2 and a a <» are then calculated. 

6. If the oc^ and a a ^ values were lower than a certain limit -got from practical experience, 

then an actuator array according to the shape of the next higher mode is created and the 
procedure was repeated from step 3. 

7. The design was stopped if the a m and values are higher that the above specified limit. 

The above mentioned limit for the control performance indices a a2 and has to be 

chosen from practical experience. As a rule of thumb values above 0.2 were considered to be 
acceptable. 

4.3 Sensors Placement 

The optimal sensor placement depends in twofold respect on the output matrix C. On 
one hand the magnitudes of the single elements indicate the sensitivity of a single sensors, and 
on the other hand the condition of C indicates their correlation. 

4.3.1 Sensitivity of Single Sensor 

The first objective is to obtain preferably high elements of C. This can be performed 
by the vivid placement of the sensors on the antinodes of the controlled modes as these are the 
points of maximum curvature. However, this approach is limited to simple structures and a 
small number of modes. In more complicated cases a sensitivity performance index for each 
sensor arrangement has to be calculated. Hence, this is carried out in the same lines as the 
control performance indices derived in (4.9) to (4.19). The only difference being that the 
placement indices now indicates the controllability of the sensor with respect to the n th 
mode. The B„ has now to be inserted as a unity matrix, since it is assumed that the actuators 

are placed ideally. The sensitivity indices a S2 , a s „, p j2 and for the sensors are determined 
like their respective actuator performance indices. 

4.3.2 Condition of The Output Matrix 

The solid computational solubility of the equation system determined by C depends 
primarily on its condition. As the same kind of sensors and measurement conditions is used 
for each sensor to obtain the velocity, the C matrix becomes C v matrix which is quadratic. 
The mathematical most precise way to solve a system of linear equations is given when its 
matrix is orthogonal. In case of linearly dependent rows of the matrix considered, the solution 
gets singular. If C v is near to singularity, the calculation of the desired modes will be 
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5. The control performance indices a a2 and oc a =o are then calculated. 

6. If the a a2 and a a <» values were lower than a certain limit -got from practical experience, 

then an actuator array according to the shape of the next higher mode is created and the 
procedure was repeated from step 3. 

7. The design was stopped if the a a2 and ot a „ values are higher that the above specified limit. 
The above mentioned limit for the control performance indices a a2 and a a ~ has to be 

chosen from practical experience. As a rale of thumb values above 0.2 were considered to be 
acceptable. 

4.3 Sensors Placement 

The optimal sensor placement depends in twofold respect on the output matrix C. On 
one hand the magnitudes of the single elements indicate the sensitivity of a single sensors, and 
on the other hand the condition of C indicates their correlation. 

4.3.1 Sensitivity of Single Sensor 

The first objective is to obtain preferably high elements of C. This can be performed 
by the vivid placement of the sensors on the antinodes of the controlled modes as these are the 
points of maximum curvature. However, this approach is limited to simple structures and a 
small number of modes. In more complicated cases a sensitivity performance index for each 
sensor arrangement has to be calculated. Hence, this is carried out in the same lines as the 
control performance indices derived in (4.9) to (4.19). The only difference being that the 
placement indices now indicates the controllability of the sensor with respect to the n th 
mode. The B„ has now to be inserted as a unity matrix, since it is assumed that the actuators 

are placed ideally. The sensitivity indices a S2 , a^, $ s2 and (3^ for the sensors are determined 
like their respective actuator performance indices. 

4.3.2 Condition of The Output Matrix 

The solid computational solubility of the equation system determined by C depends 
primarily on its condition. As the same kind of sensors and measurement conditions is used 
for each sensor to obtain the velocity, the C matrix becomes C v matrix which is quadratic. 
The mathematical most precise way to solve a system of linear equations is given when its 
matrix is orthogonal. In case of linearly dependent rows of the matrix considered, the solution 
gets singular. If C v is near to singularity, the calculation of the desired modes will be 
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impossible because this will lead to mathematically unprecise results. Therefore a scalar 
quantity to classify the status of C v between orthogonality and singularity is needed. This 
scalar quantity can be found by computing the condition number of the matrix considered. 

The condition number is calculated as 

cond([CJ) = , (4.20) 

W min([Q]) 


where w max , w min indicates the highest and lowest singular value of the matrix C v , respectively. 
The singular values are computed using a Singular Value Decomposition (Numerical Recipes 
[20]) An orthogonal matrix has a condition number of one, a singular matrix has a condition 
number of infinity and a arbitrary matrix a real number in between. 


A measurability index can thus be formulated which includes the condition number and one 
of the developed sensitivity indices a or p. The measurability index of the sensor array iy is 

. _ mag , (4 

'' cond([C ]) 

where mag s is the magnitude index a or p. The index ranges between one for the best 
possible configuration and zero for a defective one. 


4.3.3 Design procedure 

The design is done by the following steps: 

1 . The mode shapes of the structure is to be found out. 

2. The modes to be controlled (n c ) are then identified. The number of these modes equals 
the number of sensors. 

3. The sensors are placed near the bending antinodes of the mode shapes. 

4. The output matrix C v and the the measurability index i s are then calculated. 

5. If the index i s is low, then the condition number and the magnitude index has to be 
checked. 

6. If one of the p indices is too low then it means one of the sensors is placed awkward. The 
line of the output matrix with the lowest values indicates which sensor is then changed. 

7. If the condition number is too high, then the column with the lowest values has to be 
found and one of the sensors must be moved to one of the antinodes of this referring 
mode. If it is to possible to identify one particular line or row which contains only low 
values then one of the sensor is moved arbitrary. 
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8. A new placement combination is then used and the steps from 4 of this procedure is 
repeated until the measurability index fits the requirements. 
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Active Vibration Control of a Beam - Simnlation and Experimental Results 


This chapter introduces the simulation and experimental results for the IMSC applied 
to a vibrating simply-supported beam. The optimum actuator and sensor placements are found 
out for controlling the modal amplitudes of the first four modes. The control problem is 
simulated for different vibration cases and the effects of spillover are inspected. The 
experimental work is carried out and the results are discussed. 

5.2 Material Properties 

The material and geometric properties of the steel beam specimen is listed in the table 

below. 


Property 

Parameter 

Value 

Unit 

Length 

1 

1.126 

m 

Breadth 

b 

0.04 

m 

Height 

h 

0.005 

m 

E-Modulus 

E 

2.1E11 

Nm' z 

Density 

P 

7890 

kgm' 3 


Table 5.1: Properties of the beam 


To control the vibration of abeam, PZT materials are used as actuators and PVDF materials 
are used as sensors. The electrical and mechanical properties of these materials are listed in 
Table 5.2. 






“ -simulation and Experimental Results 


Property 

Parameter 

PVDF | 

PZT 

Unit 

Thickness. 

t 

110E-6 


m 

Piezo strain Constant 

d 31 

23E-12 




mv 

Electromechanical 
coupling factor 

*31 

12 

34 

% 

Y oungs-Modulus 

1 E 

3E9 

6.66E10 

MKSBSOm 

Relative permitivity 

£/£<) 

12.5 

2400 



Tabic 5.2 : Properties of PVDF and PZT materials used in the experiment 
5.4 Modal Parameters 


An experimental modal analysis is carried out in order to determine the modal parameters 
namely eigen frequencies (to), mode shapes(vjr) and the modal damping ratios(Q of the 
structure under consideration. The values determined from the Experimental Modal Analysis 
(EMA) closely resemble the computed values calculated from the analytical closed form 
solution. The table below gives the modal information of the first four modes desired to be 
controlled. 


Mode 

No: 

Eigen 

frequency [Hz] 

Damping ratio 

[%] 

Eigen frequency 

(analytical) [Hz] 

1 . 

10.27 

11.43 

9.2256 

2. 

37.65 

3.86 

36.902 

3. 

83.49 

1.4 

83.0304 

4. 

149.49 

1.18 

147.609 


Table 5.3 : Modal parameters of the beam specimen 


5.3 Actuator Placements 

The PZT actuators are placed based on the actuator placement indices and control 
performance indices. If the values of a fl2 and a a ^ are higher than a satisfactory limit, in this 

study the limits are set to a a2 >0.4 and a a „ >0.2, the positions of the actuators are accepted. 
Further, the frequency response functions (FRF) of the beam, excited by a shaker, are 
compared with the FRFs obtained by exciting the beam with the actuators. Thus this is a good 
method to check the actuator positions selected by using a good control performance index. 
Table 5.5 shows the values of various control performance indices for the selected actuator 
position, which are shown in fig 5.1. 
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PZT actuators 



Fig 5.1: Placement of Actuators on the beam 
The table below gives the coordinates of the actuators placed on the beam. 


Actuator Nr. 

Xj Coordinate [m] 

x 2 Coordinate [m] 

1 

0.133 

0.203 

2 

0.529 

0.599 

3 

0.819 

0.889 

4 

0.903 

0.973 


Table 5.4 : Coordinates of the actuators 

The control performance indices for this actuator location was obtained as: 


^Ol 

0.41600 

(X a oo 

0.26559 


0.28440 

Pg°° 

0.18575 


Table 5.5 : Control performance indices for actuator locations 
The figure below gives the comparison of the Shaker excited and actuator excited FRF for 
sinusoidal excitation. The program to calculate the FRF is based on (2.39). 
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5.4 Sensor Placements 

The placement of PVDF sensors are based on both the calculation of the sensitivity 
indices and the condition of the output matrix [Chapter 4], The initial sensor locations for 
iteration were selected to be the points of antinodes of the control modes. The positions 
were changed to get the best condition and hence the best measurability index . The table 
below gives the values for the sensitivity indices, condition and the measurability index for 
the selected sensor location shown in fig 5.3. 


PVDF point sensors 



Big 5.3 : Sensor placement on the beam 
The table 5.6 below gives the coordinates of the sensor positions on the beam. 


Sensor Nr. 

Xj - Coordinate [m] 

x 2 - Coordinate [m] 

1 

0.5625 

0.5775 

2 

0.3525 

0.3675 

3 

0.2525 

0.2675 

4 

0.8325 

0.8475 


Table 5.6 : Coordinates of the sensors on the beam 
The performance indices and the condition number for the selected sensor placements is 
outlined in table 5.7 below. 


a a2 | 0.4369 


0.32808 

P- 

0.3419 

p.- 

0.2383 

Cond [C] 

8.7085 

h 

0.0501 


Table 5.7 : Performance indices of the selected sensor placements on the beam 
Photographs of the beam specimen with actuators and sensors bonded to it is shown in 
Appendix 5. 
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5.5 Simulation 

A program was written down in Matlab in order the simulate the modal filtered 
feedback control of a simply-supported beam for both the transient and the steady state 
excitation. The time-modal displacement plots of the beam with and without control and the 
effects of spillover are simulated. The Runge-Kutta fourth order integration, which is 
outlined in Appendix 2, is used for time integration. 


5.5.1 Transient vibration 

The free vibration of the beam with and without control, after displacing the beam in 
the middle by 5 mm is studied. The figures below gives the time-modal displacements plot 
of the first and fourth modes with and without control. For control, the time lag for 
application of control is also taken into consideration into simulation. 

4 4 

3 

2 

1 

modal amplitude 

I q I [m] 0 

[a] 

-2 

-3 

0 100 200 300 400 500 600 700 


X 10 



modal amplitude 
I q I [m] 


[b] 


time [1/1000 s] 




Without control 
With control 



_ 2 I i i i i i i 1 1 1 

0 50 100 150 200 250 300 350 400 450 



time [1/1000 s] 


► 


Fig 5.4 : Simulated Modal displacement vs time plot of [a] first and [b] fourth mode 
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The largest contribution to the response w(x,t) of the beam comes from the lowest 
mode. Control spillover does exist. By including the effects of fifth, sixth and seventh 
modes into the previous algorithm, the spillover effect is able to be observed although it is 
very small. The figure below gives the time-modal displacement plot of the seventh mode 
when the vibration of the structure is controlled by focusing control of the first four eigen 
modes. The increase in modal displacement q is a clear indication, that this particular mode is 
excited at at t-20 ms as a result of the controller. 


modal amplitude 
I q 1 (m) 



time (1/1000 s) 

Fig 5.5 : Spill-over effect on the seventh mode 


5.5.2 Steady State Vibration 

The vibration of a simply-supported beam by exciting it by a harmonic excitation of 
magnitude IN is studied. The figure below gives the modal displacement vs time plot of the 
fourth mode when the structure is excited at the fourth eigen frequency. The excitation is 
removed after the vibration reaches steady state condition and the die down of vibration is 
simulated with and without control. 



Fig 5.6 : Simulated modal displacement vs time plot for fourth mode [a] without and [b] 
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5.6 Bonding Procedure 

The bonding of the PZT and PVDF layers on the beam is one of the most important 
detail which has a large impact on the performance of the piezoelectric as an active damping 
material. In this section a brief note on the bonding of actuators and sensors have been 
touched upon. 

5.6.1 PZT Actuators 

The PZT materials are obtained in a standard rectangular size 75x25 mm. After 
marking the positions of their placement, they are bonded to the structure taking an epoxy 
glue. This glue has a good resistance to shear forces and provides good metal - PZT bonding. 
The beam PZT contact is established by means of silver ink glue which provides a good 
electrical contact so as to ground the structure. The anti-symmetric wafer configuration, see 
chapter 2, is used for the experiments. The poling directions are adjusted as shown in fig 
5.5. Here a voltage is applied to one electrode of the PVDF-actuator, where the other 
electrode is grounded. 



E = Electrical field 

fig 5.6 : Bonding of PZT actuators 


5.6.2 PVDF Sensors 

PVDF point sensors of size 1.5 x 1.5 cm are bonded to the structure as shown in the 
figure below. The sensors are covered by an aluminum foil and are connected with shielded 
cables to the signal amplifier in order to shield the signal from interfering with the 
surrounding electro-magnetic field. The outputs of the sensor are then connected to the 
amplifier which then combines both the output signals and amplifies them. 
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Double sided sticky band for bonding 

Fig 5.7 : Bonding of PVDF sensors 

5.7 Experimental Set-Up 

The experimental set up for the vibration control of the beam is shown in the Fig 5.7. 
The set up consists of two independent circuits. One circuit is used for the vibration control 
of the beam , whereas the other circuit is used to monitor the resulting vibration. 



Cables 

Fig 5.8 : Experimental set up for vibration control of the beam 
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The main part of the controller circuit is a four channel digital signal processor (DSP). 
After the control algorithm is generated on the PC and downloaded to the DSP, the DSP 
works stand alone with a sampling frequency of 20 kHz. The input channels of the DSP are 
connected via sensor amplifiers to four PVDF sensors. The PVDF material is mono oriented 
with e 31 = 4pC/N and e 32 = 0 . The amplification factor of the sensor amplifier is set to about 

100. The controller board converts the analog signals into digital ones, calculates the modal 
velocities, integrates the modal velocities, multiplies the state vector with the gain matrix, 
calculates the retransformation by using the inverse input matrix and converts the digital 
result into analog signals. These signals are amplified by four high voltage amplifiers. The 
maximum amplifier voltage is limited to 100 volts so as to be well below the breakdown 
voltage of PZT, which would otherwise damage it. The output of all high voltage amplifiers 
is connected to the four PZT actuators. 

The monitor circuit is independent from the control circuit. An electrodynamic shaker 
is used as vibration source and is integrated in the monitor circuit. ...The excitation force is 
sensed by a piezoelectric force transducer at the mounting of the gear box. The response of 
the structure is measured with two accelerometers at two different locations of the beam. 
These locations are selected to obtain significant peaks in the FRFs at the four controlled 
modes. The acceleration response form the accelerometers is obtained by writing a program 
on a HP-VEE machine which helps in generating sinusoidal excitation, sweeping through the 
frequency range 0-150 Hz and to record the acceleration response at each frequency. 

5.8 Experimental Results 

The control of the beam vibration is carried out After inserting the right gain 
coefficients. A good reduction of the vibration amplitude was observed after control. The fig 
5.8 shows the modal velocity vs time plot with / without control for the fourth control mode at 
148 Hz. The table 5.7 below shows the amount of reduction obtained in the amplitude of 
vibration for each mode as obtained from fig 5.9. 


Mode No. 

Frequency [Hz] 

Reduction [%] 

1 

10 

40% 

2 

37 

62% 

3 

83 

28% 

4 

149 

74% 
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time [s] pr- 

Fig 5.8 : Modal velocity vs time plot for the fourth mode 

The above plot shows that there is a good reduction in the modal amplitudes for the 
fourth mode at 149 Hz. The acceleration vs frequency plot for the beam with and without 
control obtained after measuring the acceleration for each frequency at a point of maximum 
displacement is shown in the fig 5.9 below. It is clear from the figure that a good reduction 
in vibration is obtained for the first, second and fourth modes respectively. The shaker 
excitation was positioned at 23 cm along the x axis. The excitation was of 7 N magnitude. 
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Chapter 6 

Active Structural Acoustic Control of a Car Structure 


In numerous industrial applications, structurally radiated noise is a persistent problem 
which is often poorly alleviated by passive means, particularly at low frequencies. In ASAC 
the control inputs are applied directly to the structure in order to reduce or change the 
vibration distribution with the objective of reducing the overall sound radiation. Though it is 
obvious that completely reducing the overall structural response with active vibration control 
would lead to the attenuation of the sound radiation, various modes of vibrations have 
differing radiation efficiencies and some modes are better coupled to the radiation field than 
others. Hence, only selected modes need to be controlled, rather than the whole response. 
This chapter gives the details of experiment and results obtained after carrying out ASAC on a 
car hood. . Fig 6. 1 shows a generic arrangement of a distributed elastic system excited by an 
oscillating disturbance. 

In the experiment the car hood is the structural specimen. The car hood can be 
considered as a plate, although its shape is quite complex. The excitation is done with a point 
load excitation applied by a shaker. The modal information is obtained after performing an 
experimental modal analysis (EMA). The better radiating modes are identified by measuring 
the sound pressure. Effort to control these modes are done with the help of PVDF materials, 
applying the theory of plates with piezoelectric materials and controlling the vibration with 
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6.1 Sound Radiation of a Vibrating Surface - The Rayleigh Integral 


A very brief introduction of the sound radiation of vibrating surfaces which is 
understood to perform the experiments is outlined here. The evaluation of the Rayleigh 
integral is probably the conceptually simplest approach for calculation of the sound field 
radiated by a vibrating flat surface which is imbedded into a rigid infinite plane. The Rayleigh 
integral gives the complex valued pressure at a given field point p( r) in terms of the complied 
velocities v(r s ) field onto the structure S. 

It is given as follows 


p(r) = 

s 2 nR 


( 6 . 1 ) 


where R = lr-r 5 l and v(r s ) is the component of the complex velocity normal to the surface S, 
while Pq is the density of the acoustic medium. Essentially the integral evaluates the sum of 
the fields of a distribution of elemental sources, each having a complex volume velocity 
v(r s )dS. Fig 6.1 shows the rectangular plate in an infinite baffle showing nodal lines, 

coordinate system (r,0,<j)) and an elemental sound source associated with the plate motion. 



Fig 6.2 : An elemental sound source associated with plate radiation 
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The radiation efficiency (o mn ) of different mode types in the low frequency range when the 
acoustic wavelength is much greater than the plate dimensions 
ka{ wave number)«l , a = b can be given as (Fuller [1]) 


for m,n both odd, 

32(ka) 2 


for m odd , n even, 

8 (ka) 4 

® mn ~ ~ 2 2 5 ’ 

3m n K 

for m,n both even, 

2 (kaf 

mn 1 5 m 2 nV ' 


( 6 . 2 ) 


(6.3) 


(6.4) 


The dependence of the radiation efficiency of these mode classes on increasingly high powers 
of ka shows that the three classes exhibit radiation efficiencies which are respectively 
characteristic of monopole, dipole and quadrupole type sources as shown in the fig 6.2 below 
which illustrates the corner monopole model of low frequency plate radiation. 


i u*. 

+ m odd, n odd + 



(dipole) 


m even, n odd 
(dipole) 


+ 


+ 


m odd, n odd 
(quadrupole) 


Fig 6.3 : Corner monopole model of low frequency plate radiation 


Thus, it is clear that the odd-odd modes contribute more to sound radiation as compared to 
the other classes. In even-even mode case there is a perfect cancellation of the neighboring 
cells thus reducing the sound radiation. In the experiment with the car hood, the modes 
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acoustic wavelength is much greater than the plate dimensions 
ka ( wave number)«l , a = b can be given as (Fuller [1]) 


for m,n both odd, 

32 (ka) 2 


for m odd , n even, 

8(faz) 4 

"" 3 m 2 rcV ’ 
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2 (kaf 

15 m 2 n 2 7V 5 
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(6.4) 


The dependence of the radiation efficiency of these mode classes on increasingly high powers 
of ka shows that the three classes exhibit radiation efficiencies which are respectively 
characteristic of monopole, dipole and quadrupole type sources as shown in the fig 6.2 below 
which illustrates the corner monopole model of low frequency plate radiation. 
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Fig 6.3 : Corner monopole model of low frequency plate radiation 


Thus, it is clear that the odd-odd modes contribute more to sound radiation as compared to 
the other classes. In even-even mode case there is a perfect cancellation of the neighboring 
cells thus reducing the sound radiation. In the experiment with the car hood, the modes 
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which contribute to the sound radiation are identified after carrying out a sound pressure 
measurement, see also section 6.5.3. 

6.2 Experimental Modal Analysis 

An experimental modal analysis is decided to be carried out for finding the modal 
system parameter, namely the eigen frequencies co„ , modal damping ratios Q and the modal 

shape- function V F„. This is the best procedure to obtain results from the car structure as it is 

quite difficult to model the geometry and the boundary conditions of the structure for a FEM 
analysis. In this sections the important concepts understood in order to carry out EMA are 
highlighted. 


6.2.1 Used Hardware 

Basically, there are three main measurement mechanisms for EMA, namely: 

1 . Excitation mechanism 

2. Sensing mechanism 

3. Data acquisition and processing mechanism 

A schematic representation of the basic hardware is shown in Fig 6.3. In this figure the 
basic set up for an experimental modal analysis is illustrated. 

The analyzer converts the analogue time domain signals generated by the transducers 
into digital frequency domain data that can be subsequently processed with the help of digital 
computers. The validity and accuracy of the experimental results basically depend on the 
understanding of signal acquisition and processing. 



FFT analyser response transducer 


□ □□□□□ 


O 
CL 


signal generator 


power amplifier 


signal conditioner 



force transducer 


exciter (shaker) 
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6.2.2 Building Blocks of the FRF Measurement 


Linear Spectrum (S x (Af), Sp(Af)) 

The linear spectrum is the frequency domain representation of the time history. 
Transformation from time to the frequency domain is accomplished with the Fast Fourier 
Transform (FFT). The linear spectrum is a complex function of frequencies with values at 
discrete frequencies determined by the frequency resolution (Af) of the measurement. The 
response spectrum Sx is given by 


SJmAf) = 


' 2 ' 


AM 

^At^x(nAt)e~ j 




(6.5) 


n=0 


Auto power Spectrums : (G^y, Gpp) 

The Auto power spectrum is the magnitude-squared of the linear spectrum and is real 
valued function of frequency. It is formed as the product of a linear spectrum S x and it’s 

complex conjugate S x * 

Gpx = S X S X * . (6.6) 

Cross power Spectrum : ( Gxp) 

The Cross power spectrum is a measure of the mutual power between two 
measurements (e.g. input force and response). It is a complex valued function of frequency 
and contains relative phase information between each measurement component. It is formed 
as the product of the linear spectrum of the response and the complex conjugate of the linear 
spectrum of the input signal. 

g XF ~ S X S F* ■ ( 6 - 7 ) 

6.2.3 Estimating the Frequency Response Function 

A traditional measurement-system model used to describe the FRF measurement when 

noise is present on both measured input and response signals is shown in the block diagram 

below. 
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Linear spectrums : 

Sf = Actual input force to the structure 

Sm = Noise, uncorrelated with the input force 

Sf = Sf+ S m - Measured force 

S x - = Actual structural response 

S x = Noise, uncorrelated with response 

S x = S x - + S x = Measured response 

Required measurement : 

It is required to measure the transfer function in the frequency domain (FRF) as 



The conventional frequency response function estimator Hj (noise on the response only) is 
determined by using the cross input-output spectrum and the input auto-spectrum. 

H, (6.9) 

Another version of the frequency response function estimator H 2 is obtained by normalizing 
the output auto-spectrum by the cross input-output spectrum: 

( 6 . 10 ) 

GfX 


As Hj and H 2 should give the same result, an indicator of the quality of the analysis can be 
defined as the ratio of these two estimators. Thus 
H, _ 

~ n rz ~ r, 

J FF '-'XX 


— j. = ° XF - G ^ = ° XF G ? f - = |2 where; 0 < y 2 < 1 

H 2 G ff G xx Grr Gyy Grr Gy 


( 6 . 11 ) 


J FF XX 


2 

y is called the ordinary coherence function. It is a normalized coefficient of correlation 
between the measured input and response signals evaluated at each frequency. It is a real 
valued function of frequency which indicates the amount of response that is linearly related to 
the input signal. Its value ranges between zero (no relationship) and one (100% relationship). 
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The cross-section of the area where the actuator is placed on the plate is shown in fig 6.6 
below. The strain distribution for the antisymmetric wafer configuration used is shown here. 



s = Cs 


Fig 6.7 : Cross-section of the plate with the actuators on top and bottom 


The actuator, which is bonded to the surface of the plate , generates the line moments m* and 
m^ along the x- and y- axis (Fuller [1]) 

m, = B p C x (H {x - x, ) - H(x - x 2 )) (H(y - y, ) - H(y - y 2 )), (6. 12) 

m^B'Cymx-xO-Hix-^MHiy-yO-Hiy-yJ). (6.13) 

The constants C x and C y represent the gradient of the actuator induced strain and are given as 


C* = 


3 h, [(/Z, + K E*Xf- + y a X-) + (h m v a -h iv E*v p K^- + v a f-)] 
W s +h iv E*) 2 -(h iu v a -h, v E*v p ) 2 ] ' 


3/z 7 [(h,; + h iv E*){-~- + V a ~) + (h iii V a -h iv E*v p )(^ + v±)] 


C. 


2[(h ii+ h h E*) 2 -(h m v a -h !v E*v p ) 2 ] 


where 


h- t = h 2 +h 2 -h 2 -h 2 , 
K = h/ + h 2 3 - h 3 -h 2 , 
K = h 3 +/z, 3 -h 2 -h 2 , 
K=h 2 -h 2 . 
and 


E* = 


E p q-Vg 2 ) 
E a ( 1-v,, 2 ) 


(6.14) 


(6.15) 


(6.16) 

(6.17) 


H(.) in (6.12) and (6.13) represents the Heaviside step function and x u yj and y 2 are the 
cartesian coordinates of the actuator boundaries. The line moments can be replaced by the 
equivalent loads as (Fuller [1]) 
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PzacM’ y) = - 


d 2 m x (x, y) d z m y (x, y) 


V 


dx 2 


ay 


(6.18) 


The modal amplitude for vibration of a plate can be written down as, see also chapter 2. 


< In = 


(0) n 2 -Q 2 )-i(Q} n 2 7]) J* 


pt A[(0) n " -Q 2 ) 2 +co n 4 rj 2 ]i l y 


J \Pzac, '¥ n (x,y)dydx. 


(6.19) 


By substituting (6. 18) in the above equation, the integral term can be written as 


Xn y 2 r- — I Xj 3*2 

J j P „ <*. y)dxdy = -B, (C„ | Wdy 


+ c. 


*1 y i 


)'i 




■ 2/t 

dy 


y 2 x 2 

( 6 . 20 ) 


-* >’! 


where, H?,, (x, >) is substituted as (x) ^(y) and B is the bending stiffness of theplate 


B 


Et„ 


p 12(1- v 2 ) 


Finally the input matrix can be created by calculating the modal input factors of each 
actuator and each mode. The input factors of an actuator are collected into a column and the 
input factors of a row are related to a mode. The input factor b n for the n th mode can be 

written down as 

x 2 . v 2 A 

j j p Met C*> y) (*> y)dydx 

b (6.21) 

n x 2 y 2 V 

* >'i 


6.3.2 Output Matrix 

The concept of point sensors allocated in an array is used in order to find the out put 
ma trix. PVDF point sensors are defined as thin piezoelectric films with small size compared 
to the wavelength of structural deformation. The modal sensor output factors can be arrived 
at after following the same procedure as outlined in chapter 2. 
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They are given by 


¥*,=- 


t. 

A '2 


ff+A/ 


-33 


Xh Vh 

V x h - v « 


)j 


x h .n 


, dydx 

3, 


( 6 . 2 : 


th 

where A is the area of the sensor and x,- p x I2 , y ;i and y, 2 are the coordinates of the i sensor. 


6.4 Interpolation of a shape-function 

As illustrated in the section 6.3 the shapefunction y/ n (x,y) is required in an analytical 
closed form in order to calculate the input and output factors. The experimental results of the 
mode-shapes are available in the form of a table. Since small non-linearity’s of the structure 
and curve fit errors lead to local deviations from the experimentally evaluated mode shapes, a 
direct use of the table is not possible. Hence an interpolation of the measured mode shapes 
with polynomials is required in order to smoothen out the local deviations and to generate a 
analytically closed form solution. 


6.4.1 Interpolation by Using Lagrangian Polynomials 

The interpolating polynomial of degree N - 1 through N points z x \ = f(xj ), z Xl =f(x 2 ), 

z Xfl =/(xyv ) is given explicitly by Lagrange’s interpolation formula 




V^ v ^ v 2 / V " v 3y 

(x, -x 2 )(x, -x 3 ) 

**V' V N J „ , ”3/ 

Z x \ + 

,.(x, -x N ) (x 2 -x { )(x 2 -Xj).... 

• • * • V' v 

....(x 2 

] (x-x l )(x-x 2 ).... 

. 


(x M -x 2 )(x N -x 2 ).. 

/ \ * 



-z x2 +- 


(6.2: 


In the above expansion, there are N terms, each a polynomial of degree N - 1 and each 
constructed to be zero at all Xj except one, at which it is constructed to be z x j . The same 


interpolation is used in the y - direction. 

xp ( )- (y-i^Xy-ys) (y-y*) z , (y-yi)(y-y 3 ) (Jzldl z + 

2n y (yi-y 2 )(yi-yi) (yi-y N ) ^-y^i-yi) (y 2 -y N ) y2 

f (y-yi)(y-y 2 ) z 

Ow-iv-i) yN 


The displacement field in x- and y-direction is composed into a product of the polynomials 
^(x) and 'P 2/) (y), 




(6.2 
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The sampling points xj and yj are selected from the measured data. In practice the 
sampling points are chosen by regarding the mesh plots of each single mode shape. The lines 
of the mesh which show the maximum displacements and no local displacements are chosen to 
be the reference lines. For each coordinate x and y one reference line has to be defined. The 
sampling points are selected from the reference lines and plugged into the interpolation 
expansions (6.23) and (6.24). The necessary number of sampling points depended on the 
complexity of the mode shape considered. 

6.4.2 Determination of Coefficients of the Lagrangian Polynomials 

The determination of the input / output matrices starts with the calculation of the 
integrals involving the form function V F. To make the integration easier, a method of finding 
the coefficients of the interpolating polynomial is used. The idea is thus to write down the 
interpolating polynomial as 

z x = 'Pi n (x)= ci + c 2 x + c 3 x 2 + + CtfX NA (6.26) 

where the c/ are required to satisfy the linear equation 

2 am 1 r i r ~\ 
x { x, .... x, c r z xl 

2 N- 1 

x 2 x 2 ... x 2 c 2 z x2 

. ‘ = . (6.27 

2 AM r _ 

% *** _ N _ _xN _ 

(This is the so called Vandermonde matrix which has the problem of being quite ill- 
conditioned.) 

Similarly it could be written that, 

'P in 0)= c i * + ci*y + c 3*/ + + c N*y N ~ l ■ ( 6 - 28 ) 

The coefficients c t and c,-* can be found out by making use of the function value interpolation 
written down in (6.23) and (6.24). If (6.23) is interpolated to find the value of the 
interpolating polynomial at x = 0, then this value is evidently q. Now q is subtracted from 

z x f and each is divided by corresponding x f . Throwing out one point ( one with the smallest 
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xi being the best candidate) the procedure couid be repeated to find c 2 and so on (Numerical 
Recipes [20]). 

After getting the coefficients and plugging in (6.25), the values of v F„(x,y) is then 

substituted into (6.21) and (6.22). The integration is then carried out and the values of input 
and output matrices can be found. 

6.5 Experiments and Results 

The experiments carried out to control the vibration and sound radiation of the car 
hood is outlined in this section. 

6.5.1 Set Up for Experimental Modal Analysis of the Car Structure 

The first and foremost part to perform the EMA on the car structure is to form the grid 
on the hood to approximate the geometry of it. The size of the hood is 130 x 138 cm. 400 
grid points are set out on the hood, with a resolution of each grid being 5x5 cm approximately. 
The geometry is subsequently build up in the post-processing software with the approximation 
that the car hood closely represents a plate structure. The structure is checked up for linearity 
before going on with the EMA. This is done using the reciprocal theorem according to which 
the FRF of the response at one point due to the excitation at a second point must be same as 
the FRF of the response of second point due to excitation of the first. The impulse hammer is 
used as the exciter. Three acceleration transducers are used for obtaining the response. 

The experimental set up is illustrated in fig 6.8. After choosing a good excitation point 
the response function (FRF) from each point is found by exciting this point with the hammer. 
Since the FFT analyzer is equipped with only four channels the measurement is repeated 
several times by changing the position of the acceleration transducers. The coherence function 
provides a check on how good the measurements are thereby helping in evaluating the 
influence of noise and uncorrelated signals. Results with coherence between 0.9 and 1 are the 
ones, accepted for further processing. The measured data are windowed with the help of a 
exponential window. 
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After measurement is complete the curve-fit process is used to obtain the modal 
parameters. The curve fit is done with various methods. The global and residues curve fit 
yields reasonable form-functions. 

6.5.2 Results of the Experimental Modal Analysis 

After carrying out the EMA on the car hood, and finishing the curve fit process the 
modal information is found. The table 6. 1 gives the values of natural frequencies and the 
modal damping ratios of the first 10 modes. 

The shape of the frequency response function is obtained after a modal peak search , 
which is an average of all FRFs at 400 points, is shown in fig 6.9. 




















Chapter 6 : Active Structrural and Acoustic Control of a Car Structure 


64 


The first four computed mode shapes are shown below. The mode shapes are found to 
be separate modes with no coupling after looking into the modal assurance criteria which 
gives a check of orthornormality of the eigen modes. 



Fig 6.10 : First four mode shapes of the hood 


All mode shapes in the lower frequency range of interest (100 - 400 Hz), are very complex in 
nature. In addition, the density of modes is high. Around 20 modes are present in the first 
200 Hz. Hence, a visco-elastic passive damping coating is applied onto the bottom of the 
hood, to decrease the density of the mode shapes. 

6.5.2.1 EMA results with Passive Damping Coat 

The passive damping coating is cut out from a large viscoelastic sheet and fixed under 
of the hood. By using a passive layer non-linearity is introduced into the structure. However, 
since the coherence function is still of good quality it is assumed that the non-linearity is not 
large, therefore EMA can be carried out. The results of the EMA shows that the eigen modes 
in the lower frequency range are still present and the modes in the higher frequency range are 
damped out. In addition the eigen frequencies in the lower frequency band are found to be the 
same as these without the passive treatment. The FRF of the hood with passive coating is 
shown below, it is clearly seen that the peaks have blunted out showing that more damping is 
present. The magnitude of FRF has reduced and the modes in the higher frequency range have 
been damped out 





Fig 6.11 FRF of the passive damped hood 

6.5.3 Experimental Set Up for Sound-Pressure Measurement 

The purpose of the sound-pressure measurement is to identify those modes which make 
a maximum contribution to the sound radiation. For this an experimental set up as shown in 
the Fig. 6.7 is utilized. Seven microphones are used to find out the far-field acoustic radiation 
of the hood. The electrodynamic shaker is used for the excitation of the structure. The point 
of excitation is chosen that the shaker excitation could excite all the modes in the frequency 
range of interest (100 - 400 Hz). 
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The more accurate sound power measurement is not carried out here. Hence the 
microphones are placed in such a manner that maximum local information could be obtain 
from various parts of the hood. Microphones are placed above the center, at the outer edges, 
in the front and back of the hood by placing them such that the information from the far field 
is obtained. The readings from the microphones helps in the identification of the modes which 
couple effectively to the acoustic media. A HP-VEE measurement system is used for 
generation of sinusoidal signals for excitation of the hood. The signal is swept through the 
frequency range of interest with a step size of 0.5 Hz. A Pi-control algorithm is written to 
ensure constant force excitation into the structure. A stand alone time of 30 seconds is set 
before the measurements from the microphones are recorded. This is done to achieve steady 
state condition for the vibration. 

6.5.4 Results of Sound Pressure Measurement 

The sound pressure measurement is carried out with and without the passive damping 
coating applied on the car hood. The graph below shows the sound pressure distribution over 
the frequency range of interest for a constant steady state sinusoidal excitation of 45 N, which 
is the average of the readings obtained from each microphone. 



frequency [Hz] i 

Fig 6.13 : Bar chart of the sound pressure measurement 
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Five modes which contribute significantly to the sound radiation characteristics are identified 
and selected for control. These five modes in the frequency range of interest are listed in 
table 6.2 . 


No: 

1 

2 

3 

4 

5 

Eigenfrequencies of 

selected modes [Hz] 

107 Hz 

119 Hz 

144 Hz 

200 Hz 

237 Hz 


Table 6.2 Modes selected for control of sound radiation 

Another interesting observation is that the effectiveness of the passive damping 
treatment on controlling sound radiation is more predominant in the higher frequency range 
and less effective in the lower frequency range. The figure below gives an illustration of 
comparison of the sound pressure measurement, which are averaged of all microphones 
conducted on the hood with and without passive coating. 



Fig 6.14 : Sound pressure with and without passive damping coat 
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6.5.5 Interpolated Mode Shapes of Controled Modes 

Interpolation of the eigen function (\j s) of all five identified modes for control is carried 
out as stated in section 6.5. Thus, the eigen function could be computed into an analytical 
form for finding the B and C matrices. The original and interpolated mode shapes are shown 
in fig 6.12. The interpolated mode shapes are more or less free from the local faults in the 
original modes. The interpolated mode shapes are accepted to resemble the original shapes 
based on the number of points in phase and error in displacement of each point. The best 
shape which shows a low value for both the aforesaid comparison parameters are used as the 
interpolated eigen function for further analytical work. The fig 6.14 gives both the original 
and interpolated mode shapes for all five modes selected for control. 

MODE 1 (109 Hz) 


measured 



interpolated 



MODE 2 (119 Hz) 


measured 



interpolated 



MODE 3 (144 Hz) 


measured 


interpolated 
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measured 


interpolated 



Fig 6.15 : Measured and interpolated mode shapes 


6.6 Actuator and Sensor Placements 

The sensor and actuator placements are decided by going through the procedure 
outlined in Chapter 4. In this section the placement of five actuators and sensors on the car 
hood is illustrated. 

6.6.1 Actuator Placement 

The PVDF actuators are placed based on the actuator placement indices and control 
performance indices. If the values of a 0 2 and a^are higher than a satisfactory limit, the 

positions of the actuator are accepted for control. Further, the FRF’s of the structure excited 
by a shaker is compared with that of the structure excited by the actuators. Thus the actuator 
positions selected by getting a good control performance index is further checked for good 
performance in controlling the selected modes. Big sheets of PVDF materials, which measure 
24x20 cm, are used to obtain larger line moments induced into the structure. Table 6.3 and fig. 
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6.15 give the positions of the selected PVDF actuators on the hood. The values of various 
control performance indices for the selected actuator position are shown in the Table 6.4. 



Xj coordinate 

x 2 coordinate 

Yi coordinate 

y 2 coordinate 

unit 

1 

-94 

-67.5 

2.3 

21.3 

cm 

2 

-94 

-67.5 

23 

42 

cm 

3 

-94 

-67.5 

50.5 

69.5 

cm 

4 

-51 

-24.5 

-15.5 

3.5 

cm 


-44 

-17.5 

29.5 


cm 


Table 6.3 : Positions of PVDF actuators on the car hood 



Fig 6.16 : Positions of PVDF actuators on the car hood 
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^a2 

0.35060 

0 C a ~ 

0.26350 

Pa2 

0.1911 A 


0.14918 


Table 6.4 : Control performance indices for actuator locations 
The above obtained values of control performance indices for the actuator position is 
the best obtained after iteration of different positions of the actuators. The fig. 6.16 below 
gives the comparison of a shaker excited FRF and actuator excited FRF for sinusoidal 
excitation. The program to calculate the FRF is based on (6. 19). The actuator positions are 
adjusted so that the shape of the FRF of the actuator position excites the modes in the same 
way as excited by the shaker. 


modal amplitude 
Iql [m] 



Fig 6.17 : FRF shapes of the shaker and actuator induced sinusoidal excitation 


6.6.2 Sensor Placement 

The placement of PVDF sensors which are bioriented ( 031=632 = 4pC/N) is based on 
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initial sensor locations for iteration are selected to be the points of antinodes of the control 
modes as these are the points of maximum curvature. The positions are changed to achieve 
the best condition and hence the best measurability index. The table 6.6 gives the values for 
the performance indices, condition and the measurability index for the selected sensor 
locations shown in fig 6.17 and documented in table 6.5. 


No. 

Xj coordinate 

x 2 coordinate 

yi coordinate 

y 2 coordinate 

unit 

1 

-75.75 

-74.25 

46.25 

47.75 

cm 

2 

-35.75 

-34.25 

-0.75 

0.75 

cm 

3 

-30.75 

-29.25 

46.25 

47.75 

cm 

4 

-20.75 

-19.25 

22.75 

24.25 

cm 

5 

-15.75 

-14.25 

46.25 

47.75 

cm 


Table 6.5 : Position of PVDF sensors on the car hood 



Fig 6.18 : Sensor positions on the Car hood 
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^a2 

0.41003 

d a co 

0.35703 

Pa 2 

0.36954 

P«- 

0.32533 

Cond [C] 

7.5725 

h 

0.54144 


Table 6.6: Performance indices for the selected sensor combination 


The photographs of the actuators and sensors placed on the car hood are shown in fig 
6.18 and in the Appendix 5. [a] shows the placement of the PVDF actuators and sensors on 
top and [b] shows the actuators on the bottom of the hood. The shaker and the point of 
excitation can also be seen in the second picture. In the performed experiments to control the 
modes of vibration, only three modes could be controlled as two of the sensors are placed in 
wrong position due to measurement mistakes. The three sensors and actuators are 
reconfigured to control the last three modes. 



Fig 6.19 : PVDF actuators / sensors placed on [a] top and [b] bottom of the hood 
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6.7 Experimental Set Up for Vibration Control 

The experimental set up to control the vibration of the car structure is depicted in fig 
6.20. As with the beam experiment the set up consists of two independent circuits, namely 
the controller circuit and the monitor circuit. 


Sensor Antplifir 



The only difference between the experiment with the car structure an the beam is that PVDF 
materials are used as actuators instead of PZT. The amplification factor of the signal amplifier 
is set to 100. The maximum actuator voltage is limited to 230 volts so as to be well below the 
breakdown voltage of the PVDF material. 

In the experiment, the number of controlled modes are reduced from five to three, as 
two of the sensors are placed wrongly during application of the sensors to the hood. Due to 
ttmo ohnrtorrA fhrpp mnHpc namp.lv the. ones at 144 Hz. 200 Hz and 237 Hz are controlled 
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after reconfiguring the input and output matrices. Actuators 1, 2, 3 and 4 are utilized for 
control of the aforementioned modes .Actuators 3 and 4 are working together to control the 
same mode. The sensors 1 , 2 and 3 are used for observing these modes of vibrations. 

In the monitor circuit, the electrodynamic shaker is used as vibration source. The 
excitation signal is sinusoidal covering a frequency band form 140 to 300 Hz. A photograph of 
the experimental set up is shown in fig 6.20 below. 



Fig 6.21 : Photograph of the experimental set up for control of Car hood vibration 

6.8 Experimental Results of vibration control 

The sensors 1,2 and 3 are found to give a good condition number of 4.5 to observe these 
three modes. The actuators 1, 2 and 3&4 are used for control of the structure. 

The results of acceleration measurements of two accelerometers placed at two 
different positions which does not lie on the nodal lines of the three modes controlled are 
shown in fig. 6.16. These accelerometer plots are obtained by exciting the structure with 
sinusoidal excitation and sweeping through the frequency range with a PI controlled force of 
4N. It is observed from the figures that the acceleration of the controlled modes are 
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significantly reduced. It is clear from the second accelerometer plot that a good reduction of 

acceleration is obtained for mode 2 (200 Hz ) and mode 3 (237 Hz), respectively. 

A 


accelaration 

[ms -2 ] 




frequency [Hz] 
Accelerometer 1 


► 

Accelerometer 2 


Fig 6.22 : Measured acceleration at two points on the hood 


A sound pressure measurement with three microphones located at the center and the 
edges of the hood is done both with and without control. The results for the three 
microphones are summed up in fig 6.24. The readings are obtained after exciting the 
structure with a controlled force of 4 N which is implemented with the help of a Pi-algorithm. 
The force variation over the frequency range is shown in fig 6.23. It is clear from the figure 
that all microphones which give the local sound pressure at their locations recorded an 
average decrease of 3 dB over the whole frequency range. Though this is a small value, a 
reduction of 5-6 dB is obtained for the controlled modes. 
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The fig 6.25 shows the modal velocity vs time plot for the control of mode with eigen 
frequency 237 Hz. This mode is the 3 rd mode for control and hence the third modal velocity 
has the maximum amplitude, thus showing that the modal filter is working properly. The plot 
shows the modal velocity - time plot both before and after control. 



Time [s] 


Fig 6.25 : Modal velocity - time plot for the third mode of control - 237 Hz 

A comparison plot of the 3 rd modal velocity of the third control mode - 237 Hz, before 
and after control is shown in fig 6.23. It is clear from the figure that there has been a 
reduction in the velocities thus showing that the structural vibration of this mode has been 
reduced. 



Time [sj 
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Effect of passive coat 

The active modal control was also carried out after applying the passive damping coat 
on to the bottomof the hood. Fig 6.27 shows the Sound pressure- frequency plot for the 
excitation of the hood with passive layer. A constant sinusoidal excitation of magnitude 9 N 
was applied on to the structure. 



Frequency [Hz] 


Fig 6.27 : Sound pressure measurement of hood vibration control with passive coat 
The accelaration - frequency plot of the vibration level is also obtained by placing an 
accelerometer on one of the most vibrating parts of the hood. This is shown in fig 6.26. 


A 


acceleration 

[ms' 2 ] 



frequency [Hz] 


► 


Fig 6.28 : Measured accelaration of hood vibration control with passive coat 
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Chapter 7 
Summary 


Modal filtered feedback has been applied to control the structural and acoustic 
response of the vibration of a car hood using discrete PVDF actuators and sensors. For this a 
pretest has been conducted by controlling the vibration of a simply supported beam. After 
modelling the beam with actuators and sensors, the modal superposition has been utilised to 
obtain the modal input, B and output, C matrices. Independent modal space control (IMSC) is 
utilised to control the first four modes of the beam vibration. 

The success of the vibration control of the beam has been extended to achieve the 
structural and acoustic control of a car hood. The modal information, namely the modal form 
functions, eigen frequencies and the damping ratios has been obtained after an experimental 
modal analysis. The modes which radiate more sound have been identified with the help of a 
sound pressure measurement set up. The mode shapes have been converted into an analytical 
closed form by doing an interpolation. The actuators and sensors have been placed on the 
structure based on the sensor and actuator performance indices. The vibration of three modes 
of the car hood in the lower frequency band have been controlled. 

A hybrid active-passive control has been investigated by attaching a viscoelastic 
coating on to the bottom of the car hood. The experimental results show that this approach 
holds much promise in reducing the sound radiation over a wider frequency range. 

The following have been achieved : 

• control of the first four modes of the vibration of a simply supported beam. 

• structural and acoustic control of three modes of a car hood. 

• active - passive hybrid control of the vibration and sound radiation of a car hood. 
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Chapter 8 : 

Conclusion 


The experiments with the car hood and the beam have shown that the modal controller 
works with a high efficiency and that vibration can be reduced successfully over a wide 
frequency band. The concept used has its limitations in the sense that the deviations by the 
curve fitting and truncation of residual modes, which reduce the accuracy and lead to sub 
optimal vibration reduction. In addition it has been observed that a hybrid passive and active 
treatment for vibration and noise reduction has larger implications in the control of sound 
radiation of a structure as it enables the sound reduction over a wide frequency band. 


Scope for future work 

It is suggested that the above work can be extended to control more than three modes 
of the car hood. This can be achieved by designing the actuator and sensor locations to control 
more modes. A better performance in vibration reduction could then be observed as more 
energy could be applied into the structure. In addition, a higher level of adaptation could be 
achieved by going in for unconstrained algorithms such as the adaptive feedforward control 
where the input to the structure evolves with time as a function of changing environmental or 
operational conditions. It would also be worthwhile to investigate the effectiveness of modal 
control over other direct output feedback algorithms such as the Positive position feedback. 
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§§ Appendix 1: 


Linear Piezoelectric Theory 

Linear piezoelectric theory couples the electric and mechanical fields. In this appendix 
the linear piezoelectric theory is reviewed. This theory is based on the quasi-static assumption 
in which the electric field is balanced with the elastic field so that these two fields can be 
decoupled at a given time instant. The direct and converse piezoelectric effects are written as : 

{D} = [e]{S} + [£ s ]{£}, (A.1.1) 

{T} = [c E ]{S}-[e] l {E}. (A. 1.2) 

where {T} is the stress tensor (second order); {5} is the strain tensor (second order); [E] is the 

E> th 

electric field vector; {D} is the electric displacement vector; [c ] is the elasticity tensor (4 
order) evaluated at constant electric field; [e] is the tensor of piezoelectric stress coefficients 

(3 rd order); [e]' is the transposed tensor of [<?]; [e 5 ] is the dielectric tensor (2 nd order) evaluated 
at constant strain. If the piezoelectric tensor [e] is set to zero, the above two equations become 
the conventional dielectric and Hook’s equations. 

There are a number of ways to write down the elastic, piezoelectric and dielectric 
governing equation, of which the strain and electric displacement is most common and 
outlined below: 

{ S}=[s E ]{T} + [d]'{E }, (Al.3) 

{D} = [d]{T} + [£ T ]{E}. (A1.4) 

where [ d ] is the piezoelectric strain constant matrix and [5 ] is the compliance matrix 

evaluated at constant electric field. 

The various relationships between the elastic, piezoelectric and dielectric constants are 


fc £ ] = [5 £ r‘ 

[e] = [d)[c E ] 

[ e s ] = [e T ]-[e][d}‘ 


outlined below: 


(A1.5) 
(A 1.6) 
(A1.7) 



§§ Appendix 2: 

Fourth Order Runge-Kutta Time Integration For Beam vibration Simulation 

In Modal Coordinates 


For simulation of adaptive structures without control feedback loop, it can be advantageous 
to use direct time integration schemes to solve the second-order equations of motion, 
Examples are Houbolt, Wilsons-0 and Newmark schemes. The existence of control feedback 
loops, especially with actuator, sensor or observer dynamics makes the direct time integration 
schemes difficult and hence Implicit and Explicit schemes based on the first-order state 
differential equation are preferred. The well know among them is the old workhorse - The 
Runge-Kutta method. 

The steps for the time integration of the Modal state-space equation is outlined below: 

• Assigning the value of the load p z on the structure. 

• Evaluation of the matrices A, B, C ,G and force vector h. 

® Determination of the initial values of the state space vector x. 

• Carrying out the time integration of equation (A 2. 1) with definite time step t. 

x- Ax + B (-G C“'x ) + hp z . (A 2. 1) 

The fourth order Runge-Kutta method requires four evaluations of the right hand side of the 

above equation per step size t. For the integration from the interval t/. to fy+j the evaluations 

can be done as shown below 
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Now the state-space vector for the (k+1) time step can be written down as : 


** + i = ** + 


Ax lk + Ax 2k + Ax 3j( . + Ax 4) , 


The integration can be continued till the required final time limit. 


(A 2.2) 
(A2.3) 

(A 2.4) 

(A 2.5) 


(A2.6) 
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§§ Appendix 3: 

Condition Of A Matrix By Singular Value Decomposition (SVD) 

Any m x n matrix A can be factored into a set of three matrices [U], [w] & [V] with special 
characteristics. The matrices [U] and [V] are orthogonal and [w] is diagonal. This factoring is 
called SVD. Thus, 

[A] = [U] [w][V) T (A 3.1) 

If the matrix A is square, N x N then, U, V, and w are all square matrices of the same size. 

Their inverses are also trivial to compute. It follows immediately that the inverse of A is 

A ~ 1 = V . [diag (1 /\Vj)\.U T (A 3.2) 

The only thing that can go wrong with this construction is for one of the wfs to be zero, or 

(numerically) for it to be so small that its value is dominated by round off error and therefore 
unknowable. If more that none if the wfs have this problem, then the matrix is even more 

singular. 

Formally, the condition number of a matrix is defined as the ratio of the largest (in 
magnitude) of the wf s to the smallest of the wf s. A matrix is singular if its condition number 
is infinite and it is ill-conditioned if its condition number is too large, that is, if its reciprocal 
approaches the machine’s floating-point precision, For e.g. less than 10' 6 for single precision 

or 10' “ for double . 
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§§ Appendix 5: 


Photographs of the beam and car hood specimens 



[a] and [b] Beam with Sensors and actuators 



[c] Piezoelectric actuator and sensor on the car hood 
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